Flow transitions in two-dimensional foams 
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For sufficiently slow rates of strain, flowing foam can exhibit inhomogeneous flows. The nature of 
these flows is an area of active study in both two-dimensional model foams and three dimensional 
foam. Recent work in three-dimensional foam has identified three distinct regimes of flow [S. Rodts, 
J. C. Baudez, and P. Coussot, Europhys. Lett. 69, 636 (2005)]. Two of these regimes are identified 
with continuum behavior (full flow and shear-banding) , and the third regime is identified as a discrete 
regime exhibiting extreme localization. In this paper, the discrete regime is studied in more detail 
using a model two dimensional foam: a bubble raft. We characterize the behavior of the bubble 
raft subjected to a constant rate of strain as a function of time, system size, and applied rate of 
strain. We observe localized flow that is consistent with the coexistence of a power-law fluid with 
rigid body rotation. As a function of applied rate of strain, there is a transition from a continuum 
description of the flow to discrete flow when the thickness of the flow region is approximately 10 
bubbles. This occurs at an applied rotation rate of approximately 0.07 s . 

PACS numbers: 83.80.Iz,83.10.Gr,64.70.Pf 



I. INTRODUCTION 



Nonlinear viscoclastic materials are observed to exhibit 
inhomogeneous flows. One source of inhomogeneous flow 
is geometric. When subjected to certain boundary con- 
ditions, an inhomogeneous stress will be applied to the 
fluid, resulting in an inhomogeneous flow. A classic ex- 
ample is a yield stress fluid in a Couette geometry (flow 
between concentric cylinders). A yield stress fluid only 
flows if the stress is above a critical stress, and in a Cou- 
ette geometry, the stress decreases as a function of the 
radial distance from the inner cylinder. Therefore, it is 
possible to generate flow close to the inner cylinder where 
the stress is above the yield stress, and at the same time, 
there will be no flow beyond a critical radius. The crit- 
ical radius corresponds to the point at which the stress 
has decreased below the yield stress [1]. Though inho- 
mogeneous flows are not a new phenomenon, it has been 
challenging to directly measure them due to the inherent 
opacity of most complex fluids, such as foams, granular 
matter, colloids, emulsions, pastes, etc.. Direct measure- 
ments of such flows has focused on two-dimensional sys- 
tems [2, 3] (for which the entire material is observable) 
and non-optical techniques in three dimensional systems, 
such as MRI techniques [4] . A surprising element of these 
direct studies of velocity profiles is the observation of in- 
homogeneous flow that can not be understood in terms 
of the geometry and the resulting stress distribution. 

Direct measurement of velocity profiles has covered a 
range of diverse systems. These include worm-like mi- 
celles [5, 6], lyotropic lamellar phases [7, 8], granular mat- 
ter [9-11], slurries and pastes [4, 12], and foams [2, 3, 13- 
16]. In most of these studies, inhomogeneous flows take 
the form of flow localization or shear localization. This 
refers to the fact that the system divides into two spatial 
regions: a flowing region and a stationary (or solid-like) 
region. The breadth of systems exhibiting this general 
behavior naturally leads to questions of universality of 



the observed flow localization. In this regard, the pro- 
posed jamming transition [17] provides a natural context 
for considering flow localization. Essentially, the jam- 
ming transition refers to the transition from solid to fluid- 
like behavior (or vice- versa) in a system as a function of 
density, temperature, or externally applied stress. Of 
particular interest is the correspondence between transi- 
tions due to a critical stress and the more familiar transi- 
tions (such as the glass transition) as a function of tem- 
perature and density [18]. If one is to consider shear 
localization in this context, it is natural to divide the be- 
havior into two categories: continuous and discontinuous 
transitions. 

The division into continuous and discontinuous transi- 
tions is based on the expected analogy to phase transi- 
tions. In the case of jamming transitions as a function of 
applied stress, the relevant variable is the rate of strain. 
In continuous transitions, the rate of strain is continuous 
as the system makes the transition from flow to no-flow. 
For this class of behavior, the velocity as a function of 
position in the system is often well-fit by an exponential 
[2, 10]. Such continuous behavior was the first class of 
inhomogeneous flows that was identified as potentially 
universal across different systems, having been observed 
in both granular and foam systems. Though not central 
to the work in this paper, it is worth noting that in gran- 
ular matter, the velocity profile has been attributed to a 
spatial variation in the density that develops during flow 
[10]. This results in a spatial variation of the viscosity 
that is the source of the shear localization. For foam, 
exponential profiles were observed in foam confined in a 
Hele-Shaw cell [2]. However, experiments [19] and simu- 
lation [20] suggest that this type of continuous transition 
is caused by drag due to the confining plates. 

More recently, experiments have identified examples of 
discontinuous transitions. Such flows have been observed 
in various slurries and pastes [4, 12], including three- 
dimensional foam [16], in two-dimensional model foams 
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known as bubble rafts [3], and in worm-like micelles [5] 
and lyotropic lamellar systems [7]. In these flows, the 
rate of strain is discontinuous across the system, often 
at a transition from shear flow to rigid body type be- 
havior. Discontinuous flows have not yet been observed 
in simulations, though a number of simulations exhibit 
continuous flow localization [20-24]. 

In the context of a jamming transition, a discontin- 
uous transition represents the natural analog to a first 
order phase transition. Therefore, discontinuous transi- 
tions represent an important class of transitions if one 
is to fully understand the implications of the jamming 
paradigm. Furthermore, within the context of rhcological 
models, these discontinuous flows require non-standard 
constitutive relations. Essentially all the standard stress- 
strain relations used to describe complex fluids assume a 
continuous dependence of the stress on the rate of strain, 
and hence, the rate of strain must be continuous by defi- 
nition for any transition from fluid to solid-like behavior. 

The above discussion assumes that the materials are 
well described by a continuum model, either in the con- 
text of a specific constitutive relation or the jamming 
transition. For the case of flow in a three-dimensional 
foam, there has been a detailed study of the flow behav- 
ior that combines standard rheological measures of stress 
with direct measurement of velocity profiles using a Cou- 
ette geometry [16]. For these studies, a large system size 
was used that allowed for the categorization of the flow 
into two classes. When the width of the flowing region is 
above a critical value, a single constitutive relation that 
is based on the existence of a critical rate of strain is used 
to describe the flow of the foam over a wide range of rates 
of strain. This is to be expected in the continuum limit. 
When the width of the flowing region is less than a crit- 
ical value, a different type of behavior is observed. This 
has been called the "discrete flow" regime. 

The discrete regime was observed to occur when the 
flow is localized to a spatial width that is less than ap- 
proximately 25 bubble diameters. This corresponded to a 
critical rotation rate for the driving cylinder (in this case 
the inner cylinder) on the order of 0.3 s _1 [16]. For the 
discrete regime, there does not exist a single continuum 
model that describes all of the flow curves. For example, 
when the system is driven by rotating the inner cylinder 
with a constant rotation rate, the torque as a function of 
the applied rotation rate does not follow a well-defined 
curve. 

Previous measurements on a bubble raft suggest that 
the discrete regime occurs in two dimensional foam as 
well [3] . A bubble raft consists of a single layer of bubbles 
floating on the water surface [25, 26]. Studies of the 
flow behavior of a bubble raft using a Couette geometry 
demonstrated the coexistence of flowing and non-flowing 
regions with a rate of strain discontinuity between the 
two regions [3]. In these studies, constant rotation of 
the outer cylinder was used to generate the flow. Two 
different external rotation rates were studied, and in both 
cases, the velocity profile in the flowing region was well 



described by a power-law fluid. However, different power- 
law fluid models were required for each case. Because 
the total system size was only 25 - 30 bubbles across, 
this is consistent with the expectation of a discrete flow 
regime. The possible connection between these flows and 
the observed discrete regime in three dimensional fluids 
provided the motivation for the study reported on in this 
paper. 

In this paper, we report on a more complete study of 
the flow in the small system size bubble raft to elucidate 
both the nature of the solid to fluid transition (i.e. is it 
continuous or discontinuous) and to determine if there is 
a transition from the continuum limit to discrete flow. 
We focused on low rates of strain and considered two dif- 
ferent system sizes. We use a foam confined between two 
cylinders, and apply a constant rate of rotation of the 
outer cylinder. We confirm the existence of a transition 
between a discrete flow regime and a continuum limit as 
a function of the rotation rate of the outer cylinder. In 
both cases, the flow can be described as the coexistence 
of a power-law type fluid and a rigid body, with a dis- 
continuity in the rate of strain. However, in the discrete 
regime, one can not use a single, consistent power-law 
model for all rotation rates. The rest of the paper is or- 
ganized as follows. Section II describes the experimental 
setup and techniques. Section III discusses the models 
used to describe the data, and Sec. IV presents the re- 
sults. Section V is a discussion of the results. 



II. EXPERIMENTAL METHODS 

The experimental system consisted of a standard bub- 
ble raft [26] in a Couette geometry (two concentric cylin- 
ders) . The system was driven by rotating the outer cylin- 
der at a constant angular speed. Stress measurements 
were made using the inner cylinder, which was suspended 
on a torsion pendulum. The studies used two values of 
the outer radius (R): 7 cm and 9 cm. The inner cylinder 
had a fixed radius of r» = 2.85 cm. The bubble raft was 
produced by flowing regulated nitrogen gas through a hy- 
podermic needle into a homogeneous solution of 80% by 
volume deionized water, 15% by volume glycerine, and 
5.0% by volume Miracle Bubbles (Imperial Toy Corp.). 
The bubble size was dependent on the nitrogen flow rate, 
which was varied using a needle valve. A random distri- 
bution of bubble sizes was used, with an average radius of 
1.7 mm and 15-22 bubbles across, depending on the outer 
radius. In previous versions of this system, the bubbles 
were generated separately and transferred into the appa- 
ratus. For these experiments, the bubbles were generated 
directly in the Couette apparatus. For each setting of 
the outer radius, essentially the same set of bubbles were 
used. Occasionally, near the end of a run, some bubbles 
would pop. Data was only used up to the point the first 
bubbles were observed to pop, and the bubble raft was 
filled in before the next run. For additional details of the 
apparatus, see Ref. [27]. 
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As mentioned, the system was driven by rotating the 
outer cylinder at a constant angular velocity, Q. The 
range of angular speeds used was 0.01 s _1 < fl < 
0.35 s _1 . The first layer of bubbles at either bound- 
ary was not observed to slip relative to the boundary. 
At the outer boundary, this was due to the curvature of 
the boundary. At the inner boundary the first layer of 
bubbles was held in place with metal fins attached to the 
boundary. Due to the finite size of the bubbles, this re- 
sulted in an effective inner radius on the order of 3.1 to 
3.2 cm, depending on the details of the system. 

As the experiments focus on the average velocity pro- 
file and the corresponding rate of strain (or shear rate) as 
a function of radial position, it is useful to review what 
is expected for Newtonian fluids in a Couette geometry. 
Due to the cylindrical geometry, the shear rate is not uni- 
form across the system and is given by j(r) = r ^^r^- 
Here vg is the azimuthal velocity of the bubbles. Because 
we rotate the outer cylinder, it is useful to normalize the 
velocity by the expected velocity for rigid body rotation 
V(r) — Or, where r is the radial position of interest. 
Therefore, we will often refer to the scaled azimuthal 
velocity v(r) — vg(r)/ilr. It should be noted that be- 
cause we rotate the outer cylinder, the bubbles in the 
outer portion always move. Therefore, with our setup, 
what distinguishes the "flowing" region from the "solid" 
region is the type of motion. The solid region exhibits 
rigid-body rotation, and has v(r) = 1. This is in contrast 
to many experiments in which the inner cylinder drives 
the system and "solid" behavior corresponds to a zero 
velocity for bubbles in the outer portion of the system. 

To measure the velocity, roughly one third of the 
trough was digitally recorded using a frame grabber. The 
time interval between images was selected so that the 
fastest moving bubbles could be accurately tracked from 
frame to frame. The radial coordinate is divided into bins 
of width 0.15 cm and 0.18 cm for the 7 cm and 9 cm data, 
respectively. This represents roughly an average bubble 
radius per bin. The details of the velocity measurements 
are given in Ref [3]. 

The other bubble motion of interest is the nonlinear 
bubble rearrangements, referred to as Tl events. A Tl 
event is a neighbor switching event that involves four 
bubbles. Two bubbles that are initially neighbors lose 
contact, and two bubbles that were not neighbors be- 
come neighbors. In contrast to the automated velocity 
measurements, Tl events were measured by stepping the 
digitized images one frame at a time and visually search- 
ing for the location and time at which Tl events occurred. 
Due to the associated motions of the other neighboring 
bubbles, Tl events are relatively easy to detect by hand, 
and for all reports of Tl events, two researchers indepen- 
dently tracked the events. It should be noted that the ve- 
locity can be determined automatically even if not every 
bubble is tracked. However, to reliably track Tl events 
requires monitoring every bubble. This is the reason for 
the need for identifying Tl events by hand. Currently, we 
are working on improving this apparatus so tracking Tl 



events is automated (as it has been done in other systems 
in our lab [19]). 

III. MODELS 

The data reported in Sec. IV will be discussed in 
the context of two standard continuum models for non- 
Newtonian fluids: Herschcl-Bulklcy model and power law 
fluid model [28] . The choice of these models allows us to 
test both the continuity of the transition between the 
fluid and solid behavior and the transition from contin- 
uum to discrete behavior. By construction, the Herschel- 
Bulkley model is continuous in the rate of strain, and the 
power law model allows for a rate of strain discontinuity. 
We will show that the power law model, with a rate of 
strain discontinuity, is the most consistent with the data. 
Once this is established, the power law model provides an 
effective method of characterizing the transition from the 
continuum limit to the discrete flow regime. In this sec- 
tion, we review the main elements of each model that are 
used to analyze the data. The details of the derivations 
are left to the Appendix, as indicated. 

The key element of the Herschcl-Bulkley model is the 
yield stress, r . If the stress is below the yield stress, 
the material acts as a solid. For stresses above the yield 
stress, the material obeys the following constitutive rela- 
tion for the average stress a as a function the of rate of 
strain 7: 

a( 7 ) = r + M7 ™ (1) 

[28]. As one can see, this model is explicitly continuous 
in the rate of strain. Also, it is important to point out 
that it provides a good fit to previous measurements of 
stress versus rate of strain in bubble rafts. Some typical 
values for a bubble raft are r = 0.8 ± 0.1 mN/m and 
n = 0.33 [3, 29]. Second, the model has a built in physical 
mechanism for the transition from solid-like behavior to 
flow: the yield stress. 

The second model was selected based on previous ve- 
locity profiles [3] and MRI measurements of velocity pro- 
files in three dimensional foam [16]. These experiments 
suggest that the system is best described by a power- 
law fluid coexisting with the elastic solid regime. For a 
generic power-law fluid, 

a(ff) = M7", (2) 

In this case, there is no specification of where the fluid- 
solid transition will occur. To account for this, a modified 
power-law fluid that explicitly includes a critical rate of 
strain (j c ) was introduced in Ref. [16]: 

^(7) = Ki/ic)" fOT 7 > 7c (3) 

for 7 < 7c, 7 — 0. Notice, this is fundamentally different 
from the Herschcl-Bulkley model in which r c is set by a 
critical stress. For this power-law model, r c is determined 
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by 7c, and there is a built in discontinuity in 7 at r c . 
Localized flow occurs when there is a coexistence of a 
flowing region for r < r c (where 7 > j c ) and a "solid" 
region for r > r c (where 7 = 0). 

Independent of the model, the stress in a Couette ge- 
ometry has the form a(r) = C/r 2 , where the constant 
C is determined by the boundary conditions (either at 
the inner cylinder or r c , whichever is more convenient). 
Combining this relation with the particular constitutive 
model, allows one to solve for v(r) = vg/flr. This in turn 
can be fit to the velocity data as a test for each model. 
For the Herschel-Bulklcy model, we get for r <r Cl 

where 

For r > r c , the solid body behavior gives v = 1; adjoining 
the solutions for r < r c and r > r c yields a smooth curve 
where the fit parameters are the exponent n, ri and r c . 
For the power-law case, we find 



v(r) = 



,2/n 



B. 



(5) 



Applying the same boundary conditions as before, A — 



{nr c f/ n /{r 2 l /n - rt ra ) and B = rf' 1 /{rf " - r c z/ ")- In 
this case, because of the rate of strain discontinuity, a 
smooth continuation of the power-law solution and the 
solid body curve does not exist. Therefore, the data in 
the range 0.25 < v < 0.95 are fit to Eq. 5 with n, and 
r c as fit parameters. Equivalently, r c can be computed 
as the intersection the line v = 1 {vg = Qr) with a fit to 
Eq. 5 with A, B and n as free parameters. Notice, j c is 
not a fit parameter, as it does not appear in the solution 
for the power-law model. To determine j c , one uses the 
resulting fit and computes j(r c ). 

In terms of testing the models for their applicability, 
it is important to note a few similarities and differences. 
In both models the three main fit parameters are n, ri 
and r c . The determination of r c and j-j provides a consis- 
tency check, as both (within a few percent) are relatively 
easy to determine by visual inspection of the data, inde- 
pendent of the selected model. Therefore, one does not 
expect to be able to distinguish the models on this basis. 
On the other hand, n has definite physical meaning. It 
gives the scaling of the stress (or the effective viscosity) 
as a function of rate of strain. For n < 1, the mate- 
rial acts as a shear thinning material. For n > 1, the 
material acts as a shear thickening material. All of the 
measurements of stress versus rate of strain suggest that 
n < 1, so this provides a physical test that the models 
must meet to be considered applicable. Another useful 
test of the models is the behavior of the rate of strain at 
r c . In this case, the models predict very different behav- 
ior by construction. For the Herschel-Bulklcy model the 
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rate of strain is continuous at r c , while for the power law 
model there is an explicit discontinuity. Here, because 
of the discrete nature of the data, it is difficult to deter- 
mine the continuity of the data by visual inspection. As 
we will show, other aspects of the fit clarify the nature 
of the transition. 

There is an additional test of the applicability of the 
Herschel-Bulklcy model: the behavior of r c as a function 
of the external rotation rate f2. Because the critical ra- 
dius is set by the yield stress, a(r c ) = t . From this, we 
can write a(r) = r r 2 /r 2 . Then as the rotation rate ap- 
proaches zero, the rate of strain in Eq. (1) does as well, 
so cr(rj) approaches tq and thus r c approaches r; t . 

The final question related to the choice of models is the 
use of these models in a discrete regime. By definition, 
both the Herschel-Bulkley model and the critical rate 
of strain power law model discussed here are continuum 
models. The key element of each is a single yield stress 
or critical rate of strain, respectively, that describes the 
properties of the foam in a continuum limit. It is difficult 
to modify the Herschel-Bulkley model to also describe a 
situation in which there is not a well-defined yield stress. 
However, if we consider the original general power law 
model, it is useful as an ad hoc model in the case of dis- 
crete flow. As we saw, the solution for v(r) in the case of 
a power-law fluid does not depend on j c . Therefore, one 
would expect the fitting procedure to work equally well 
whether the system is exhibiting continuum or discrete 
behavior. However, the results can be used to distin- 
guish between the two regimes based on the behavior of 
7c = 7(?"c) an d n. For a system in a well-defined contin- 
uum limit, these are expected to be material properties 
that are independent of the external rotation rate Vl. If 
one observes a strong dependence of these parameters on 
fl, this would represent a break down of the continuum 
limit. 

In summary, by considering two standard models for 
complex fluids (power-law fluid and Herschel-Bulkley) we 
will be able to test two features of the flow. First, we 
will be able to distinguish between transitions that are 
continuous or discontinuous in the rate of strain. Second, 
we will be able to test for a transition from the continuum 
limit to discrete flow. 



IV. RESULTS 

Figure 1 illustrates the fitting procedure for the two 
models for a rotation rate of 0.09 s _1 and outer radius 
R = 9 cm. The power law fit is illustrated in Fig. la, 
where the horizontal line indicates the rigid body rota- 
tion. The fit is continued past the critical radius to high- 
light the fact that this model gives a non-zero critical rate 
of strain. This represents a discontinuous transition. In 
contrast, the fit in Fig. lb is for the Herschel-Bulkley, and 
it is continuous by construction. From these fits alone, 
it is difficult to distinguish between these two models, 
especially at this rotation rate. Therefore, one needs to 
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FIG. 1: (color online) (a) A velocity profile (with r c — 5.9 cm) 
as a function of radial position for the system with outer ra- 
dius of 9 cm and rotation rate of 0.09 s _1 (black symbols), 
with power law fit where n = 0.515 (blue line), (b) Same 
velocity profile (black symbols) , with numeric fit to Herschel- 
Bulkley model where n — 1.056 (blue line). 
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consider more carefully the parameters for the different 
models, as they arc determined from the fits. 

Figure 2 provides a summary of the behavior of the 
critical radius as a function of system size and rotation 
rate. For comparison, we show both the calculation us- 
ing the numerical fit to the Herschel-Bulkley model (open 
symbols) and the fit to a power law/rigid body coexis- 
tence model (closed symbols). The results for r c based 
on each method are in reasonable agreement. This pro- 
vided an important consistency check on both methods. 
It should be noted that once r c « R, the various methods 
of determining r c break down. To indicate where this oc- 
curs, we still plot r c , but define it to be equivalent to R. 
In other words, r c — R is the condition that the entire 
sample is flowing. 

Three important features are highlighted by this figure. 
First, the critical radius is roughly linear in the applied 
rotation rate. (The straight line in the figure is a guide 
to the eye.) Second, for rotation rates less than 0.15 s _1 , 
changing the system size does not alter the location of the 
coexistence between the two regions. Once the critical 
radius is r c = 7 cm, in the R = 9 cm system, the critical 
radius continues to increase with rotation rate. Finally, 
the critical radius approaches r = 4.5 cm as fl goes to 
zero. The fact that r c — > 4.5 cm as ft — > s -1 is the 
first evidence for the breakdown of the Herschel-Bulkley 
model, as discussed in Sec. III. As a further test of the 
behavior of r c , we have considered two slower rates of 
strain: l.Ox 1CT 3 s _1 and 3. Ox 1CT 3 s~ x . Due to the finite 
lifetime of the bubbles, the results for the velocity profiles 
were noisier at these extremely slow rates of strain than 
the data studied in detail in this paper. However, these 
profiles clearly exhibited a r c that was close to but greater 
than 4.5 cm, consistent with the results reported in Fig. 2. 

As discussed, the exponent n in both the Herschel- 
Bulkley model and the power-law model has physical sig- 
nificance. Figure 3 summarizes the values of n obtained 



FIG. 2: (color online) Comparison of the measured critical 
radius as a function of rotation rate for outer radii of 7 cm 
(black squares) and 9 cm (red circles) . The open symbols cor- 
respond to calculation of r c using a Herschel-Bulkley model, 
and the solid symbols correspond to calculating r c assuming 
a power-law fluid coexisting with a solid-like state. 



in the various fits to the velocity profiles. Figure 3a is 
the results for the power-law model, with the R — 7 cm 
and R = 9 cm indicated by squares and circles, respec- 
tively. Figure 3b is the results for the Herschel-Bulkley 
model. The dashed line in each figure represents the 
median value for the exponents. There are two striking 
features in Fig. 3. First, neither fit provides a completely 
consistent value for n. However, the variation in the fits 
for the power-law model is smaller than the case of the 
Herschel-Bulkley model. Second, for the power-law fits, 
one consistently finds n < 1 , but for the Herschel-Bulkley 
model, one consistently finds n > 1 (though there are 
some isolated cases for which n < 1). 

As the combination of Figs. 2 and 3 effectively rule 
out a Herschel-Bulkley model, it is necessary to further 
probe the applicability of the power-law model. This al- 
lows us to confirm whether or not the transition is truly 
discontinuous in the rate of strain and to look a transi- 
tion from continuum to discrete flow. One approach is 
to consider the scaling of the velocity profiles. In this 
regard, there are two different scalings that are partic- 
ularly useful. These are presented in Fig. 4, and these 
represent the central measurement of both the disconti- 
nuity in the transition and the transition from contin- 
uum to discrete flow. First, it is useful to directly test 
for the existence of a consistent critical rate of strain in 
the context of the power-law model. To do this, we fol- 
low the scaling arguments presented in Ref. [16] . Because 
we are rotating the outer cylinder, we first subtract out 
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FIG. 3: (color online) (a) The exponent n obtained from fit- 
ting the velocity curves to a power-law fluid coexisting with a 
solid-like state versus the external rotation rate for R — 7 cm 
(black squares) and R = 9 cm (blue circles), (b) The expo- 
nent n obtained from fitting the velocity curves to a Herschel- 
Bulkley model versus the external rotation rate for R — 7 cm 
(black squares) and R = 9 cm (blue circles) . 



the rigid body rotation. This gives a new transformed 
velocity: v = (Qr/r c ) — vg/r c . In this form, v = at 
r = r c . Therefore, rewriting the rate of strain in the 
form Mr) = r-f ^ = ^ - ^, we see that the slope of 

i \ / dr r dr r 7 ^ 

v at r = r c is the critical rate of strain j c . Therefore, if 
there is a single 7 C for the material, then v as a function 
of r/r c will collapse onto a single curve. 

Figure 4a presents v as a function of r/r c . This illus- 
trates two features. First, for sufficiently high rates of 
strain, the data suggests a single value of j c . However, 
the scaling breaks down below a critical value of the ex- 
ternal rotation rate. Second, this form of the velocity 
highlights the discontinuity in the slope at the critical 
radius. 

Given the apparent dependence of j c on Q for slow 
rates of strain, it is useful to consider an alternative scal- 
ing of the data, given in Fig. 4b. Here we plot vg/Q,r c 
as a function of r/r c . For this case, we do not subtract 
the rigid body behavior, which is apparent as the linear 
regime for r/r c > 1. This was done to confirm the con- 
sistent rigid body behavior for all rotation rates. This 
scaling allows us to focus on the data for the slowest ro- 
tation rates, which did not scale in Fig. 4a. In this case, 



FIG. 4: (color online) (a) Plot of the scaled velocity profiles 
v — (Qr/r c ) — vg/r c versus r/r c . (b) Plot of ti#/Or c versus 
r/r c . Rates of strain are indicated in the legend. 



we observe a collapse of the data for slow rates of strain, 
indicating a dependence of 7 C on ft. 

The final measurement of the discrete/continuum flow 
transition is the direct measurement of 7 C as a function 
of J7. This is plotted in Fig. 5. In this case, we fit the 
data for v near r c . Two features of are illustrated by the 
data. First, for rotation rates below £1 = 0.07 ±0.02 s _1 , 
the critical rates of strain are dependent on the exter- 
nal rotation rate, but independent of the system size. 
Interestingly, the behavior is consistent with a linear de- 
pendence on f2. This is strong evidence for the break 
down of any continuum description of the flow. Second, 
we observe a cross over to a regime in which j c is in- 
dependent of SI. This occurs for value of the external 
rotation above fl = 0.07±0.02 s _1 and is consistent with 
the scaling of v presented in Fig. 4a. From this data, we 
find 7 C = 0.07 ± 0.01 s _1 , where the error represents the 
standard deviation of the measured values of the critical 
rate of strain, for our bubble raft. 

Given the possible breakdown of the continuum ap- 
proximation, it is worth testing the short time behavior 
and determining the approach to steady state flow. This 
was done by considering a range of time intervals over 
which to compute the average bubble displacements, and 
the corresponding average velocity profiles. For short 
enough time intervals, we were able to compute a his- 
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FIG. 5: (color online) Critical rate of strain as a function 
of external rotation rate for R = 7 cm (blue triangles) and 
R = 9 cm (black squares) . 



FIG. 7: Illustration of the convergence of the measured value 
of r c as a function of the time interval over which the average 
velocity is computed is increased. 
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FIG. 6: Probability distribution of measuring a particular 
value of r c based on averaging the velocity over finite time 
intervals. Four different time intervals are illustrated for a 
rotation rate of 0.07 s _1 with R = 7 cm. The time intervals 
and mean for each distribution are indicated in the figure. 
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FIG. 8: (color online) The Tl event rate for R = 9 cm and 
Q, = 0.11 s _1 is shown as a function of r (red histogram, 
left hand axis). This is compared to 1 — v(r), where v(r) = 
Ve(r)/Qr is the scaled velocity (black points, plotted against 
the right hand axis). 



togram of the probability distribution for r c . Here, the 
probability distribution is computed as follows. For each 
independent time interval, an average velocity profile is 
computed. As reported in Ref. [30], these profiles tend 
to be highly nonlinear, but there is a well-defined radius 
at which the profile deviates from a rigid body rotation. 
This point is taken as the critical radius for that real- 
ization of the velocity profile. This is computed for each 
independent segment of data from a single run, and the 
probability distribution is generated from this set of data. 



Figure 6 presents the probability distribution for the 
case of fi = 0.07 s" 1 , and for four different time intervals 
used to compute the velocity profiles. As expected, the 
longer the time interval, the narrower the distribution. 
However, even for relatively short time intervals, the full 
width of the distribution is only on the order of 1 cm, or 
about 3 bubbles. As one increases the averaging time, the 
mean of the distribution remains constant to within 1 %. 
Another way of considering the approach to steady state 
is illustrated in Fig. 7. This shows the computed critical 
radius for different sets of data using a fixed time interval. 



8 



The plot illustrates that for time intervals greater than 
50 s, the computed value of r c does not change and that 
different measurements of r c give the same value. (Due 
to the length of the run, there are two points for the 
At = 60 s data, but one point for each of the larger 
values of At.) 

A final question of interest is the spatial distribution 
of the Tl events and the correlation between the velocity 
profile and the Tl events. This is shown in Fig. 8 for Q = 
0.11 s -1 and R = 9 cm. Plotted here is a histogram of the 
number of Tl events in each radial bin (normalized by 
the area of the bin) . This is compared to a plot of 1 — v(r) . 
The velocity is plotted in this form to better correlate the 
Tl and velocity distributions. The solid line indicates 
the position of r c . One observes that the distribution 
of Tl events is consistent with the measurements of the 
distribution of r c for short time averages. As expected, 
one observes Tl events for values of r close to, but slightly 
greater than, r c . However, there are no observed Tl 
events for r sufficiently greater than r c . 



V. DISCUSSION 

The main focus of the measurements in this paper is to 
test carefully two aspects of the transition from solid to 
fluid behavior in a slowly driven foam. First, whether or 
not the transition is continuous or discontinuous in the 
rate of strain. Second, whether or not there is a transition 
from a discrete to a continuum flow regime. We used 
two standard models of complex fluids to interpret the 
experimental results: the Herschel-Bulkley model and a 
power-law model with a critical rate of strain. First, 
we will discuss the evidence against the Herschel-Bulkley 
model. 

As discussed in Sec. Ill, the initial evidence against 
the Herschel-Bulkley model is the fact that r c does not 
approach ri as f2 goes to zero (see Fig. 2). Additional 
evidence against the model is the results for n presented 
in Fig. 3. The fact that we essentially always measure 
an n > 1, in direct contrast to the measurements of 
the stress as a function of rate of strain, rules out the 
Herschel-Bulkley model. Finally, the scaling of the ve- 
locity in Fig. 4 that strongly suggests the existence of a 
discontinuity in the rate of strain is behavior that is not 
possible in the context of the Herschel-Bulkley model. 
One might be concerned that despite this evidence, one 
appears to be able to fit the data to the Herschel-Bulkley 
model (as in Fig. 1). However, given the discrete na- 
ture of the data, the apparent fit to the Herschel-Bulkley 
model is most likely a result of the non-physical values 
of n that are obtained. 

In contrast to the Herschel-Bulkley model, the fits to 
the power-law model, and the scaling in Fig. 4a, strongly 
support the applicability of the power-law model and the 
corresponding discontinuity in the rate of strain. There- 
fore, it is reasonable to state the the localized flow in 
the bubble raft is best described by a coexistence of a 



power-law fluid and a solid region, with a discontinu- 
ity in rate of strain at the coexistence point. This has 
important implications for the jamming transition, as it 
suggests the equivalent of a first order transition. Also, 
it raises the important question of the mechanism that 
determines the critical rate of strain. This will be the 
subject of future work. One promising direction is to use 
a parallel shear cell. In this case, one expects a uniform 
stress across the system, and the global rate of strain is 
set by the boundaries. This can be used to further test 
the nature of the critical rates of strain. For example, in 
this case, it was determined to be 7 C = 0.07 ± 0.01 s _1 . 

An interesting open question is which features of the 
flow are determined by the details of the bubble raft. For 
example, how does 7 C depend on the specifics of the so- 
lution used to make the bubbles, the polydispersity, the 
bubble size, etc.. Likewise, we observed that the value 
of n has potentially significant variation. For example, 
even though n < 1 based on fits to the power-law model 
(Fig. 3), it did differ from some past measurements in 
bubble rafts [3]. This difference is not surprising given 
that the details of the bubble rafts differed to some degree 
in terms of the exact nature of the bubble size distribu- 
tion and the solutions used to make the bubbles. 

The transition from a continuum regime to a discrete 
regime is confirmed by the results presented in Figs. 4 
and 5. Based on the scaling arguments from Ref. [16], 
the data in Fig. 4a is consistent with a single continuum 
model based on a power-law fluid with a critical rate 
of strain for sufficiently high rotation rates. The break- 
down of the scaling indicates the transition to the dis- 
crete regime. Directly plotting the critical rate of strain 
for each rotation rate and system size, we find that the 
transition between discrete and continuum behavior oc- 
curs at il = 0.07 ± 0.02 s _1 . By considering two system 
sizes, we were also able to establish that the behavior is 
not dependent on system size. Conversely, we observe 
that the systems we have considered are large enough to 
exhibit behavior consistent with a continuum limit. 

It is useful to compare our measurement of the transi- 
tion to discrete flow with the reported measurements in 
a three-dimensional foam. Using the results in Fig. 2 for 
the critical radius, a transition at £1 = 0.07 ± 0.02 s _1 
corresponds to a critical thickness for the flowing regime 
of approximately 10 bubbles. For comparison, the same 
transition is observed in three-dimensional bubbles for a 
thickness of 25 bubbles and a rotation rate of 0.3 s^ 1 [16]. 
At this point, it would be useful to consider both two 
and three dimensional models that capture this disconti- 
nuity to determine the source of the differences in these 
measurements. It might be dimensionality, but given the 
small differences, it is more likely related to details of 
the bubbles, such as bubble size, polydispersity, and sur- 
factant composition. As current simulations only predict 
continuous transitions, work is first needed to elucidate 
the mechanism for the discontinuous transition. 

The existence of the discrete regime also has implica- 
tions for the jamming transition paradigm. The jamming 
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transition is proposed in terms of an externally applied 
stress. However, this stress is often applied to the sys- 
tem using a constant external rate of strain. The system 
is then assumed to pass through the jamming transition 
as the stress increases as a function of total strain. It is 
not clear how the jamming transition should be under- 
stood in the discrete regime where the fluid-solid transi- 
tion does not appear to be well-described by continuum 
models. This will be an interesting direction for future 
studies. 

Finally, it is worth commenting on the insights gained 
by consideration of the short time behavior of the system. 
We focused on measurements of r c as a function of the 
averaging time (as shown in Figs. 6 and 7). This was 
useful because it provided confirmation that our results 
represent the steady-state of the system. In addition, the 
comparison of r c and the spatial distribution of Tl events 
are very suggestive. 

The direct measurement of the distribution of Tl 
events confirmed the solid-like behavior for r > r c in 
terms of the complete absence of Tl events in the bulk 
of that region (Fig. 8). However, the Tl events do ex- 
tend beyond r c by a few bubble radii. This is consistent 
with the width of the distribution of r c values for suffi- 
ciently short time measurements (Fig. 6). As Tl events 
occur right near r c (the steady-state average value), this 
will generate short-time bubble motions that can extend 
up to a few bubble diameters beyond r c . These bubble 
motions can generate deviations from rigid body motion. 
Therefore, the fluctuations in r c have an important con- 
nection with the statistics of the Tl events at the transi- 
tion from the fluid-like to the solid-like behavior. It will 
be interesting in the future to compare this to other fluc- 
tuating transitions between regions with different types 
of flow such as those that have been observed in other 
complex fluids [6, 7]. 
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VI. APPENDIX 

In this Appendix, we present some of the details of 
the derivation of the velocity profiles. The derivation of 
the profile for the power-law fluid is standard [28]. One 
combines the the general expression for the stress in a 
Couette geometry and the relation between the stress 
and the rate of strain: 



and combining all the constants, this simplifies to 

7 = D/r 2/n . 
Using the definition of 7, one gets 



(7) 



. = dvj(rl = 2/ 
dr r 



which can be directly integrated to get (dividing by ft) 

(9) 



~ v(r) — -777- - B. 



ftr 



r 2/r, 



As given in Sec. Ill, the boundary conditions determine 
A and B in terms of n, r c , and rj. 

The derivation for the Herschel-Bulkley model follows 
along similar lines, only now it is useful to explicitly write 
out the constant in the stress relation: 

cr(r) = C/r 2 = a{n)r 2 /r 2 = T r 2 c /r 2 . (10) 

Notice, that there are two equivalent ways to write the 
constant because it is based on the equivalence of torques 
in the radial direction [28] . Again, we equate this expres- 
sion for the stress with the constitutive relation. 



This gives 



T r 2 Jr 2 =T + Lij n - 



7=(r / Ai ) 1/ "[r c 2 /r 2 -l] 1 ^ 



(11) 



(12) 



In this case, integration does not result in an analytic 
expression, instead we get 



Biff ->)'*<>+<>■ ™ 



Here the constant C — because vg/r — at r = r,. 
Requiring solid body rotation vg = ftr at r = r c gives 



ftr c 



a(r) = C/r 2 = »(j/j c ) n , 



(6) 



Therefore, converting to v(r) = v$(r)/(Qr), 

*>-ff0(7)'- 1 )""'" (15) 

where 

This is a useful form for numerically fitting the data by 
approximating the integral. 
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